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Abstract 

We consider non-linear stochastic field equations such as the KPZ 
equation for deposition and the noise driven Navier-Stokes equation 
for hydrodynamics. We focus on the Fourier transform of the time 
dependent two point field correlation, <$>k(t). We employ a Lagrangian 
method aimed at obtaining the distribution function of the possible 
histories of the system in a way that fits naturally with our previ- 



ous work on the static distribution. Our main result is a non-linear 
integro-differential equation for <&k(i), which is derived from a Peierls- 
Boltzmann type transport equation for its Fourier transform in time 
$k,w That transport equation is a natural extension of the steady 
state transport equation, we previously derived for $k(0). We find a 
new and remarkable result which applies to all the non-linear sys- 
tems studied here. The long time decay of <&k(t) is described by 
$k(£) ~ exp (— a|k|£ 7 ), where a is a constant and 7 is system depen- 
dent. 

1 Introduction 

Non-linear field equations describing a system with friction, non-linearity 
and a driving noise have received much attention. Two prime exam- 
ples are the KPZ equation [j]J for the height of a non-linearly deposited 
material, 



dh 

di 



+ vV 2 h + (Vh) 2 = f 
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and the Navier- Stokes equation 

^ + ,V 2 u + u-Vu-^ = f (1.2) 
at p 

where p is determined by V • u=0. 

The noise, /, is considered Gaussian, stochastic driven, with 

(f(v,t)f(r',t'))=H(v-r',t-t') (1.3) 
We write the equations in a general form in Fourier transform: 
-7^- + Mc + ^2 Mk i l h i hl + S ^ k J lm hjfahm-- = fk(t) , 

3,1 j,l,m 

(1.4) 

where we write h as a scalar, but a minor elaboration of notation cov- 
ers vectors, e.g. M^^u?uJ . M is taken independent of the origin, and 
therefore contains #k,j+l- Note that we are using here box normaliza- 
tion, namely the Fourier transform is defined as = -7= J h(r)e lk ' r d 3 r, 
where V is the volume of the system. Consequently our M's are order 
of 1/VV and the iV's are order of 1/V etc. If ( |1.4|) is also Fourier 
transformed in time: 

iu;h ku) + z/ k h kcJ + ^ M kujiaW h i(J hi0 H = f kw , (1.5) 



where M now also contains 5u,a+6- In earlier papers, we approached 
the steady state solution of the system ( |TTT|) by deriving the Liouville 
equation for the probability P (h,i), 



BP B 

~kT ~ E Bh7 + E -f)P = ° (1-6) 

3,1 



dt ^ dhk 

which when averaged over /, taking this to be the noise with 



H = D (r - r') 6 (t - t') -»• D 0k 5 (t - t') , (1.7) 



satisfies the well known form 



dP d I d 

¥ "E^ ^ + ^ + EW:|^0 (U 

k \ j,l 



P is now the average over / of equation flOp , and the steady state 
satisfies: 



B I B 

E dj^ ( Dok dh^ + ^ k + E M whh 1^ = ° (1-9) 



The approach to equation ( [L.9|) is to derive a transport equation 



based on a self consistent method, i.e. suppose that the system can 
be developed about the model [g, |3j : 

zM Dk dk + "* hk ) p °=° (L10) 



i.e. 



P = Mexp 



where 



(«-k) = 0k (1.12) 

(1.13) 



0k being the true two point function. In Peierls' treatment of non- 
linear crystal electricity ||], 0k appears as the number of phonons nk 
and satisfies the Boltzmann equation. In turbulence 0k is the energy 
in the mode k. In granular deposition there is no name for 0k but 
perhaps we can call it the "flucton" since it measures the surface fluc- 
tuation. Peierls could use perturbation theory to derive the kernel 
of his Boltzmann equation, but since the non-linearity dominates our 
problem we need both 0k and cu^. Pq is the approximate solution 
of equation ( |1.8|) , which starts the self consistent expansion, which is 



given in j2|, |3j symbolically by: 
P = J 1- V Mhhhlj, Y MM(Kkhh),^) hh terms 

+ E MMMm/MhWMW*) + similar terms + ...\ 

(1.14) 

pictures of theses terms are given in Appendix A below. It was shown 
in refs [@, |3j that the conditions and (|1.14j ) lead in second order 

to a Peierls-Boltzmann (P.B) [Hfl form 



, _ /" M k jiMj k -iMj ^ _ f M k jiMi k; _j0 k 0i rf g^, 



|M kj i|^j^] , 



d J j = A>k (1-15) 



J CJ k + + ^1 

where it turns out that the coefficients in the terms ^k^j and have 
the effect after integration of opposite signs, and those of 0k0j and 0k0i 
are equal. A similar equation has been derived by Bouchaud and Cates 
|5j]. Note that the value of (hkhjh\h m ) is given by (h^hj) (h\h m ) + 
(hkh\) (h m hj) + (/ik^m) (hjh\) and several terms like M 2 3 / ^u; plus 
terms like M 4 4 / (X]^) 3 so that in the present theory, the four point 
correlation is not the sum of two times two point correlations. One 
should note that within this paper we are deriving an expansion rather 



than a closure. (A systematic graphology for the higher terms is given 
in appendix A. Feynman diagrams do not give the Peierls- Boltzmann 
equations derived here.) Hence a Peierls-Boltzmann [|J structure has 
emerged: 



A remarkably simple scaling argument emerges if we continue the 
series ( |1.14j ) and derive from it a series for the correlation function 
0k = {hkh-k)- The series amounts to a systematic expansion in 
"(model - reality)". Each order in the expansion will have a leading 
power in k which is 2a — T — 2/i + d greater than the previous term 
where dimensionally M = k a , (j) = k~ r , uj = k^, in d dimensions, be- 
cause symbolically, the power series is in terms of M 2 (j)/u 2 . In order 
that all terms have the same leading power, it follows that 



A well known example arises in the Kolmogoroff dimensional analysis 
of turbulence, caused by a source near k = (call this case KNS). 
There, d = 3, a = 1 and Kolmogoroff argues that T = 11/3, so that 
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2a — V — 2^ + ^ = and fi 



2a + d-T 
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the scaling argument gives \x = 2/3 (of course Kolmogoroff invokes 
the dimensional argument to obtain 2/3, but our point is that there 
is a much more general argument relating the time scale u k 1 with the 
correlation function 0k). Note that the ease with which the scaling 
argument can be checked to all orders in our expansion confirms the 
value of the method. We can see that higher terms in the expan- 
sion cannot alter powers in ( |1.1Y| ), only front factors. We develop an 
equation for to below in (3.22), but it is to be realised that whatever 
equation is deduced, all it gives is a front factor; the power is deter- 
mined by scaling. What is needed is a transport equation that can 
naturally produce behaviour that is more general than a decaying ex- 
ponential, and to do this one must treat time, or its Fourier transform 
u, as a natural extension to four dimensions, i.e. (/ikw^k'o/) • To do 
this we study the whole history distribution V ([/&kcj])- Such functions 
are of course well known in quantum mechanics after their original 
introduction by Dirac j6j , in the form: 
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Our self consistent method works for the extensions of 0k, <^k i-e. 
^k,^, ^k,w, and is presented in the next section. Several papers are 
present on this problem in the literature [|7fl-[[T^| but our method has 
the advantage of producing simple equations (simple considering the 
complexity of the problem) which allow us to produce explicit solutions 
due to the ability to check scaling relations to all orders in a systematic 
expansion. 

2 The Lagrangian formulation: a model 

In this section, we study the simple case of a single degree of freedom 
obeying a noise driven linear equation 



The main reason for considering such a simple model is that we aim 
at obtaining a first order differential equation in time for the non-linear 
systems we consider, that will match the static equations derived in 





or in Fourier variables, 




(2.2) 
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our previous papers. This task is somewhat complicated by the fact 
that the noise in equation or in any physical system cannot be 



instantaneous, since it originates in physical processes. Consequently, 
in any physical system, the time derivative of, say, (h (0) h (t)) at time 
t = is zero, and therefore a first order differential evolution equation 
cannot evolve the system in time. To understand what is going on 
and to obtain the correct matching condition, we study the system 
(|Z.1| - [2.2|) by considering a non instantaneous noise described by the 
correlation: 

(U-.) = -y + lT 1 (2-3) 

7T 



so that 



H/ir 

\hujh-w) = 7 ^— 9 i 721 ' (2.4) 



i.e. 



H /e 



-v\t\ n -l\t\ 



<t{t) = (h{0)h(t)) = ¥ — i [— . (2.5) 



Then 



$(0) =(7^ (2 ' 6) 
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and, as expected, 



j> (0) = 0. 



However, if / >> z/, i.e. almost instantaneous noise, 



•<0) = £^ 



and 



$ (r) = -D c 



for 1// << r «\jv. This is described in figure 1 



SHOW FIGURE 1 HERE. 



In the limit / — > oo, the Fokker-Planck equation for the 



of finding h at t satisfies: 



di-dh[ Do dh + u \ 



p = o, 



which gives 



+ z/$ = * > 



and 



d<S> 
"dt 



+ z/$ = £ < 
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and $ (0) = Dq/u. The awkwardness of is removed by 

putting in the full dependence on /, but more simply, as described 
above, confining ourselves to t > 0,we have the first order differential 
equation with the initial condition $ (0) = Dq/u, that implies 

a finite derivative $ (0) = —.Do at £ = 0. This matches the static 
equation obtained from the Fokker-Planck equation ( ^.1U|) 

D -i/(hh) = (2.13) 

The form of the Fourier transform of $ (t), $ (cj), suggests the struc- 
ture of <f>k,uj in the non-linear field theory. 
For the simple linear case, 

u ~ p + i/ 2 ] [u 2 + I 2 } = (2 ' M) 
where Q = [iu + z/] [zu; + /] . 

Thus the two decays of $(£), for £ > 0, are present as zeroes of ft 
(poles in in the upper half of the complex cu plane.) In the limit 
/ — > oo, the situation is similar, but there is only a single decay 

^ = ^T^ 2 = W (2 - 15) 

12 



where Q = [iuo + v\ . 

The natural model to try for the $ of the non-linear equation is 

* ku = (2.16) 

where = D -^ ^ > with no singularities in the complex to plane, 
and where gives in the first self consistent approximation a simple 
decay described by 

Q kw = zu; + u; k . (2-17) 

However, there will be a much more complicated time dependence in 
the full $ than one simple decay. (Equations ( |2.14j ) and (|2.1b|) form a 



simple example where the decay is given to first approximation by a 
simple decay but indeed the behaviour is more complicated) 

One possible definition of Q is to use the response function and 
define Q by 

^r- = (/-k-.M (2.i8) 

and employed in mode-mode coupling studies [§, [15], |T6j. We will take 



the view that we can write $ in terms of a sum of exponential decays 
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which can be extended to cover continuous distributions i.e. branch 
cuts rather than simple poles. Confining ourselves to poles for the 
moment, if we write 

"^ = E^rV= v (2-19) 

we find that is an even function in cu (a polynomial for a finite 
sum of simple decays), and it has no singularities. 
Thus 



ttk. = n t. r (2.20) 

1 li {iv + CJ k) iJ 

and 

«L = rr 7 ^— ^ (2-21) 



When $k (t) is obtained from ( |2TTg ) only the poles cj = zcjk,z con- 



tribute to £ > and their conjugates for t < 0. 

The strategy we will adopt is to construct a transport equation 
for in which the first order approximation ( |2.1V|) for £\ w will be 
useful and it will result in a higher order approximation for and 
consequently for $k (t) that will show now a decay that is much more 
complicated than a simple exponential. 
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3 The expansion 

The starting point is equation ( |1.5|) , where in order to make our no- 
tation less cumbersome, we denote the d + 1 vector (k, u) by K and 
write the equation in the form 

n 0K h K + Y< M kji/U/*L - /k = (3.1) 

We define next P {Hk', Ik} to be the distribution of the h^s in the 
presence of a given noise, /k- P is to be averaged eventually over the 
noise. The Liouville equation is now replaced by 



P = (3.2) 



JL 

which is similar to a Fermi supplementary condition. Equivalently, 
equation ( |3.2| ) can be replaced by: 



J2 <^J (tooithK + J! M kil h 3 h L -f K ) jP = 0, (3.3) 

to obtain the correct hierarchy of field correlations. This is achieved 
by multiplying the above equation by products like /iLi---^LiV and in- 
tegrating by parts. A simple example is obtained by multiplication by 
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/iK that yields the correct average of equation ( p7T|) . 



^ok (h K ) + ^ M kjl (/ij/i L > - / K = 0. 

J,L 



(3.4) 



An alternative derivation of equation ( |3.3| ) starts with consideration 



of a d + 2 dimensional system in which the equation of the form 



§- s h K (s) = 



^ok^k(s) + ^2 M JKhhj(s)h L (s) - /k + #k(sJ 



JL 



(3.5) 



is considered, where ^k(s) is a d+ 2 dimensional noise obeying 



< 0k(») >= and < 5 , k(s)s'-k(s / ) >= T> K 5(s - s'). (3.6) 



(Note that /k does not depend on s. Therefore it plays the role of 
quenched randomness). The Fokker-Planck equation for V{1ik,s}, 
the distribution of a given configuration {hj^} is given by 



K 



5s 



(9 



^KTTT 1" ^OK^K + Mkji/lj/lL - /k 

JL 



P. (3.7) 



The 'steady state' (s independent) V in the limit where X>k is zero is 



the distribution of HkjP and equation ( |3T3|) is recovered 



In order to construct an expansion for P, we write equation ( |3.3|) 
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as: 



^ d ( Ok d n ^ \ ^ d 



+E 



K 



_d_ 

dh 



K 



(^OK — ^k) h~K 



Mkji/ij/iL - /k 

JL 



P = (3.8) 



Notice that since the sum is over all K, both Qk and Q-k will appear 



in the second derivative. 



We expect an expansion, of the average of P over the noise, around 



the Gaussian 



Pq oc exp 



h\ch 



K'i K 



2 ^ $K 
K 



(3.9) 



and, as before, associate in (|3.8| ) a notional A to the second term on the 



left hand side of equation (|3.y| ) and A 2 to the third term. Expanding 
to second order in the A "Chapman and Enskog" expansion, we get 



kjl 



M k flhjh L h„ K 



+ ... 



(3.10) 



The condition that < /^k^-k > calculated to second order in A is 
equal to $k (that is the zero order result) yields an equation for $ in 
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terms of Q. 



i(Q 0K + n ,-K)$K-2^ 



M kj iMjK,- 



JL 



E 

JL 



\M^L ^ = El 



L^K 



1 1 

+ 



n K a 



K 



(3.11; 



We can recover the structure of a transport equation, familiar from 
the static case, i.e. "un-pick" £1ok + ^o,-k and 1/^k + 1/^-k by 



returning to equation (|3.1|), multiplying it by h-K and averaging over 



the distribution (|3.1U|) . Using also equation ( gTTgj) we obtain 



M kjl M JK , 



JL 



$L$K - ^ 



|M k jlP 



JL 



$J$L = 



(3.12) 

At this point we use, just in the non linear term, the first order 
approximation Qk = ico + cok- This has the advantage that now 
^j + ^l + ^ k = cjj + cji + cjk- Fourier transforming back and recall- 
ing that the zeroes of Q-k are in the lower complex cu plane, we find 



for t > the local equation 



at 



+ VkMt) - 2 V - 



CJj + CJi + CJk 



l*k(t) 



E 



IM, 



kjll 



$j(*)$i(*)=0 (3.13) 
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and 



at cjj + cji + cj k 

lM kjl | 2 



V — ®\(t)$i(t) = Ofort < 0. (3.14) 

^cjj +cJi + cj k JW w v ' 

Notice that the above is possible because the coefficients, M k ji, do not 
depend on the fourth components of the vector. The initial conditions 
with which equation ( |3.13| ) has to be solved are 



$k(O) = k , for all k, (3.15) 

where k is the static correlations. The static equation determining 
0k is 

^ M kjl M jk; . ^ |M kjl | 2 
^k - 2 > , , ^ k ~ 2^ — — | 0j0i =ko • 3.16 

It can be shown that 

f l ™ + ^^ = " k0 ' (3 ' 17) 
is an exact relation, obeyed by the exact two point function. Using 
this general result, we see that as t — ► 0+, the evolution equation 
for $ k (t), eq. ( |3.l3|) fits exactly onto the static equation for k , eq. 
( |3.1(3|) . The equations ( |3.13| ) and ( |3.1t)|) were originally derived for 

19 



the driven Navier-Stokes equations |7j, but although equation ( p. 
can give the Kolmogoroff [[T7j spectrum with a good value for the 
front factor, the boundary condition was not understood at that time, 
and this hindered further development of this approach to non-linear 
equations by this route for several decades. Note that the simplicity 
of the basic equations (|3.13H3.l5p and (gTTgj) sets this approach apart 



from mode-mode coupling theories. We have found a plausible way 
(the structure of Q) which leads to these manageable and transparent 
equations. The amazing feature of the self-consistent approach is that 
the time-dependent equation has an explicit and local dependence in 
time. It now offers a way to complete the system of functions, 0k, 
in a satisfying way. To the best of our knowledge the direct evaluation 
of the indices of 0k and cjk and the universal structure of the time de- 
pendent correlation functions are not available from other treatments. 
The simple structure of our equations now offers a way to complete 
the system of functions, 0k, <^k, in a satisfying way. We define tu^ 



20 



customarily to be given by 



which is a natural definition, if we think about a single mode decay. 
Integrating equation (|3.13fl over time, taking into account ( |3.1Y| ) and 



the fact that J °° dt = —0k, we obtain 



^ M kj iM jk ,-i ^ |M kjl | 2 
u k = v k -2) J - — J — — -0i - > J - x 

[Uj + Ul + CJkJ ^"^ [Wj + UJ\ + CJkJ 



101 / ^k 



(3.19) 



where 



H $i(t)$}(t) dt = ^ . (3.20) 
Jo u il 

Neglecting Dofc and v k in equation ( |3.1(j| ) (the static equation) and 



( p,19|) (the to equation) we find that in the inertial range 

u k = y i^m^-Ki-^) (3 . 2r 

^ (u k + + CJl) k 



A similar equation has been derived by a different method by Edwards 
and McComb ||, who used it to derive the Kolmogoroff front factor, 
achieving a reasonable value. Details of the alternative method are 
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given in the book of McComb [jTTfl - ^ * s eas y ^° see that eq.( p.l9|) 



leads just to the scaling relation discussed in the introduction because 
lu^i scales as u. This scaling relation together with the static equation 
gives 

u k = Bik 3 / 2 in KPZ in 1+1 D 

(3.22) 

= B 2 k 17 in KPZ in 2+1 D [2] 
cj k = B 3 k 2/3 in KNS [13] 

4 A closer look at the steady state in the inertial range 

The structure of the steady state equation has the form 

J A«(j, k)^ k d d j - J A( 2 )(j, /c)0j0 k _j d d j = k0 -^ k k . (4.1) 

The kernel A^ 2 ^ stemming from the |Mkji| 2 term is positive definite 
but A {1) (j,k) may attain positive as well as negative values depending 
on the specific problem and the values of j and k. For KPZ with 
noise, that is white (in space) it was shown explicitly [10] that the 
exponent characterising the leading small k dependence is obtained 
by equating the left hand side of eq. ( |Q|) to zero. The concept of 
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inertial range is given thus a definite meaning since the exponent does 
not depend neither on the source nor the 'viscosity'. For KPZ with 
noise given by ko = D$k 2p [M it was shown that up to some threshold 



Pq the inertial range still exists and the exponent characterising the 
small k behaviour does not depend on the source or the 'viscosity' 
(and therefore, its value is the same as in the white noise case). For 
p > po? the system is driven by the noise and even the leading k 
dependence depends on the source term. The problem of turbulence 
seems to combine both behaviours in an interesting way. The simplest 
conceptual picture of the inertial range is that offered by Kolmogoroff, 
where random forces put energy into a fluid at low k, and viscosity 
removes it at high k. We can make an extreme model of this by taking 
the input to be e5 (k) and the output to be which we write 

symbolically as e5 (k — oo), so that 

J A (1 Vk0j d 3 i - J A( 2 Vk-j0j rf 3 j =e[S(k)-8(k- oo)] . (4.2) 
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The Navier Stokes M. guarantees that the integral over k of the left 
hand side vanishes for any 0k 

J (^J A (1 Vk0j d 3 j - J A( 2 Vk-j0jh d 3 ^ d 3 k = (4.3) 



and so matches 



J d 3 k e [5 (k) — 5 (k — oo)] = . (4.4) 



To clarify fi4.2[) - (|4.4|) we offer a model in Appendix B. 



We have computed the value of the left hand side of equation ( fl~2|) 
for a range of values of T, and calling the quantity Z (T), we find that 
the value of T for which Z (T) = is (see figure 2) 

T = 3.6667 (4.5) 

to the accuracy of our computation (For KPZ, we have already found 
the value 2.59). 

show figure 2 here 

We have used the 0k = k~ u ^ 3 and integrated the left hand side of 
eq. (|4.2|) over a finite k sphere, interchanging the k and j integration, 
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we obtain, as expected, a non-zero result, in spite of the fact that for 
any finite k the integrand is zero. 

This results gives one the confidence to proceed to the much more 
difficult problem of the time dependence. 

It is interesting to consider briefly another example which is Navier- 
Stokes driven white noise (call it WNS). In the problem, the viscosity 
can again only influence very large k, but now the source plays a vital 
part. Integrals converge, with the solution of 

J A (1 Vk0j - J A (2 Vk-j0j = A) (4.6) 
where Dq is now a constant, and one finds 

(j) k oc ££T 5/3 L>o /3 , and 

(4.7) 

cu k oc Fk 5/3 Dl /3 

If one tried to solve this in the KPZ style of ignoring Dq, one would 
of course get Kolmogoroff again, but Kolmogoroff 's —11/3, 2/3 regime 
will not satisfy equation (fCT|). 
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5 The time dependent correlation function 



So far we have obtained the solution for the steady state correlation 
function, (h^ (t) /i_k (£)) in the two cases of KPZ and NS. Turning 
to the time dependent function, $k at first sight we are faced with 



a formidable set of equations, given by ( |3."13|) , and flXJjj ), of 



which is central. We give firstly a crude argument and relegate 



a full discussion to Appendix C. If we write $k = (fi^e rk( ^ we can 



manipulate ( |3.13| ) and ( |4. 1| ) to obtain 



5Tk 

dt 



J a 2 (i _ e -m)-Tu-^+r^ d3j (51) 



and argue (to be justified at greater length in Appendix C) the solution 
is dominated by the that part of j space where 



1 = e - r j- r k-j+ri< 



(5.2) 



i.e. T$(t) + = I\(£). We can scale out t and have 



ij + r k _j 



r k 



(5.3) 
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where I\ = r k (l) etc.. Consider this guess: I\ ~ k 2 . Then we 
would need k 2 = (k + j) 2 + j 2 which by Pythagoras's theorem gives 
the locus of j to be a sphere of diameter k. Inside the sphere 

r k > r k j + Tj (5.4) 

Outside the sphere 

r k < r k j + Tj (5.5) 

so the guess is inadequate. The only way to find a region of j space 
wherein (|5.3|) obtains is when T is linear when j equals ak such that 



1 ^ a ^ 0. The integral goes outside that locus, but the corrections 
are logarithmic and derived in Appendix C. Returning to T k (£) 

r k (t) = a |k| t x/ ^ (5.6) 

where a is a constant of the appropriate dimensions. For fi > 1, 
we derive in Appendix C the long time decay of $ k (t), including the 
logarithmic corrections 

d 1 

$ k (f) oc (alklt 1 ^) 2 exp(-a|k|t 1 ^). (5.7) 
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6 Some speculations 

Working to a kernel of order M 2 , the "transport" equation we have 
derived has the classic form with operators A, B such that 

A<f> + B<f></> = D; (6.1) 

but the solutions for different physical cases seem radically different. 

For KPZ, it appears that A and D are not central to the inertial 
range away from k = and k = oo, and the solution is given by taking 
a power law /c~ r , and deriving in detail 

B<fxf> = Z(V) (6.2) 

so that the value of V must be such that 

Z(V) = (6.3) 

and this gives the solution in the inertial range. The solution is then 
improved by adding the effects of A(f) and D. However, for problems 
rich in dimensional parameters, e.g. instead of M^i being a combina- 
tion of powers, it could be a function like e~( k +( - k+ ^ +] ) c , containing 
a new constant c with dimensions (length), likewise for A and then 
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equation is more of a standard Boltzmann form and, for example, 



a perturbation approach would be 



A- B AA + - (6 ' 4) 



which can be improved in various ways. Navier-Stokes seems to have 
aspects of both of these, for at first sight, if one tries a power, then 
B(f)(f) = has the Kolmogoroff solution derived just as the KPZ is 
derived. In fact the lack of uniform convergence puts Kolmogoroff 
into the form: 

B<fxf> = D (6.5) 

and the A(p term only comes in to polish the solution at high k. For 
WNS an attempt to use flop inevitably gives the Kolmogoroff indices 



which are quite incorrect for this problem. One must use ( |b.5|) since 
D is a constant for all k, not just being non-zero at k = as in KNS. 
Suppose now that we proceed to the next orders, which will give 

A(J) + B(J)(f) + C(J)(f)(f) + E(j)(f)(f)(j) H = D (6.6) 
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Can we now try <j) ~ k r and obtain 

z B (r) + z c (r) + z E (r) + • . • = o. (6.7) 

For KPZ, Zb (r) = gives a T which cannot be obtained from 
dimensional analysis. More detail is given in Appendix D. There is no 
reason to believe that if we proceeded to equation ( |(j.7| ) it will not give 



again a T, not quite the same as that of equation QET21), but if the self 



consistent approach works well it will be close. Certainly numerical 
simulations equivalent to all orders in (|tj.tj|) do agree remarkably well 



with the second order, (|tj. '4 ). However, in other equations, in partic- 
ular KNS, the deficiency in dimensions of the equations leads to the 
Kolmogoroff solution both by looking at Zb (r) = 0, and by balancing 
the source term with the energy cascade. What happens if we go to 
the higher orders in KNS? Could it be that equation ( 15771 ) has a non 
Kolmogoroff solution? The calculation of C, E etc. in equation ( |Q|) 



is a formidable undertaking and, unless some new dimensional quan- 
tity is needed to make C, E convergent, there seems to no reason that 
Kolmogoroff should not again be the solution. 
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8 Appendices 



In quantum field theory the algebraic series which, in for example the 
basic papers of Schwinger |nj, can be rather impenetrable, becomes 



much simpler with the use of Feynman diagrams. In this paper, we 
are proposing that the Boltzmann equation is the correct target for 
equations like KPZ, we should offer a graphology to simplify the 
appreciation of the formal expansion. This has in fact been done in 
the original study of NS many years ago |7], [T7j], and extended to 
many body problems by Sherrington [|2D| , and then fully extended 
to turbulence problems QZTj. However, it is not well known and so 



we reproduce it here, and extend the analysis to some new cases, in 
particular the expressions for fourth moments. We use the notation 
of the steady state problem, although it works equally for the time 
dependent case. The problem is to find P which satisfies f|1.8|) . P is 
to be expanded about Pq, which satisfies ( |Q|) , and this is effected by 
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introducing D and u, so that writing 

AD = D-D (A.l) 
and Au = u — v (A. 2) 

— [D— + uh-AD—-Auh + Mhh) (Po + ^i + iV--) = 
Oh \ Oh Oh J 

(A.3) 

and ascribing "A" to Mhh and "A 2 " to AD and Acj, we expand 



P = P + Pi + P 2 + --- (A.4) 

p i = G iL Mhhp o 

Oh 

p 2 = G—MhhG—MhhPo 
oh Oh 

( d d d \ 



where G is the Green function defined in ( |A.tj|) below. To evaluate Pi, 



P2, ■ • • we are repeatedly faced with the problem of finding J, where 
J Gk abc ..h' a KK ...Poll dti = J (h) (A.5) 
or since symbolically 

^M D M +uh ) G=m (A - 6) 
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and defining G by 

d ( d \ 
Khhh ...P = — (D— + uh) JP 

= P °{ D l- h ) {-l) J = P ° dj ^ 
we will show that the significant problem (significant to order V~ x rel- 
ative to other terms), is when K a fj C Ji a hbh c . . . has none of a, 6, c, . . . 
paired i.e. 6, c • • • ^ —a. Then we try J = J a bch a hbh c ... we find a 
solution provided 

Jabc... = u a + u b + u c + . . . (A. 8) 



by direct substitution. The second derivative always gives zero, and 

d 

UaK^7-h a = u a h a (A.9) 
dh a 

If ever we do find h\h-i we replace it by 4>\ + (h\h-i — 4>\) and the 
bracketed term, as in all field theories leads only to terms of order 
V' 1 . Thus as the series develops, a term in h adds an h, a-^ removes 
an h, a G inserts a term like u for each h. The first correction to 
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Pq, P\, illustrates the process: 



P 1 = ^GM k5l h i h l --Po 



k 



[GM^hh-M Po(A.lO) 
Y^PoGMwfHhh-k/fa (A.ll) 




Suppose that we wanted the value of (h a hbh c ), then it is given to 
order M by 

kji_j_j — z _^ permutations! . (A. 13) 

w k + cjj + cV ' Jl 

(Note that M has £kj+i) At this point we introduce the diagrams. 
Draw these pictures: 




J 



d 



Mhjhi 



(A.14) 



dh k 



1 
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_d d__ 



Da = 



d d 
dh^ dh-, 



-AD 




(A.15) 



_d_ 
dhx 



-Z/ k /i k = 



d 




(A.16) 



The rippled line will always point in the direction to the left, for in 
final evaluations integration by parts gives non-zero values only when 
the J^- finds an h to its left. The full lines representing h can point 
left or right. Note that cu^ = u;_k and that as in normal field theories, 
strictly speaking arrows need to be attached to the lines. 



-k 
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Our problem now is: 




Upon integration, Pq joins k and — k to give the line 



— = 0k (A.18) 

k -k 

Thus our series for P is 




(A.19) 

The non linear coupling Mkji used so far depend on the size of the 
system. In fact Myi = -^mkji where is of order 1. In the integrals 
appearing in the following we redefine the M's to be identical to the 
m's. Namely, the M's that appear in the following are of order 1. 

The condition 0k = (h^h-k) as expressed in eq. ( |1.15|) becomes 
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+ + 




(A.20) 



to which we make this commentary: 0k = 0k by definition, comes 
from < Mhhh > for odd averages must be zero, or alternatively there 
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represents 



lM kj il0iM jk) _i 3 . 
; k / ; ; a J, 



(A.2i; 



CJk + CJj + CJl 

(In the integrals we write for simplicity of presentation /, meaning 
/ = k — j as shown in the diagram. This is the convention we adopt 
in the following) where 1 (unity) is the result of 
when integrated by parts, and (J^to) 1 comes from the Green function 



i.e. h-§- h 



as in ( |A.8| ). A useful way to find this factor is to put a vertical line 



at each point where a Green function originally lay and put an to for 
each line crossed including these lines that yield unity: For example 
in the diagram above, 



M 



KM 



we obtain a factor of (cjk + cjj + Ui) . 
Likewise, 





j 






k-j 


~i 


< — J 



M 
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represents 



+ + w\ 



(A.22) 



The representation (|A.2U|) gives the steady transport equation ( |1.15| ). 



We now proceed to give a brief account of higher orders. Firstly we 
can continue ( |l.lo|) to the next order, which must be M 4 as M 3 will 
give zero. There are as many terms, so we just give typical terms. The 



four M's give 




G G 




and two typical terms are 




MhjiMj mn M nj -_ m Mjk,-i0i0m0k 



(cjj + U\ + CJk) {u>l + Wm + Wn + ^k) (^j + ^1 + ^ 



d 3 ld 3 m 



(A.23) 
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MkfiMj mn M n j- m M-k-j-i(f)i(f) I n(f)j 



and cross terms with AD and Acj e.g. 



d 3 ld 3 m 




(A.24) 



M kj iM_ Jb ,_ j ,. z 0i0j(AD) j ^ 
(cjj + cji + cjk)2cjj 



(A.25) 



Next we consider the value of (h^h-jh-i) which has the first ap- 



proximation as above in (|A.13|) . 




+ permutations 



1 M kj i0j0i 



permutations 



(A.26) 



Corrections appear at order M 3 typically 
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f ^■kmnM—j — m M i— n 

010m0n ,3 
!y J (iUj + tUi + CJ k ) (u\ + CJ n + CJ ) (cjj + U\ + CJ m + CJ n ) 

(A.27) 

Finally the four h correlation is given by 




(hbh-jh-ih-m) = (Mfy^j-m + <fat<l>}8k,m5j,-l + <M\$kj&l t -m + 



1 



-m,—n, 



-k<Pk<Pl<Pn 



+ 



V (Wk + CJj + CJi + CJ m ) (Wk + w m + CJ n ) 

J_ M_ J ^ ;n M n _ fc;m k 0i0 m 

V (cJk + CJj + CJl + CJ m ) (cJk + w m + w n ) 



(A.28) 
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and many similar terms, multiplying vastly at the next M 4 order. 
Notice that the series can be thought of as an expansion in 

M 2 <t> d 

and as has been noted earlier, a scaling relationship essential for the 
plausibility of the expansion, requires this combination to have a be- 
haviour like k°. For the time dependent equations the graphology 
holds with the 3 vector k replaced by a four vector, but with one 
other change which is the interference of the external noise with the 
non linear terms, however this does not have an effect until a higher 
order is reached than the second order we have used in this paper, and 
so we will not discuss it further. 
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B 



At first sight it is unusual for an integral operator to give 5 functions, 
but this can be illustrated by a simple model. Consider: 



r dy [y__ f°° dy fE = Y(x) fB n 

Clearly 

/•OO 

/ Y (x) dx = (B.2) 
Jo 

for one needs only to interchange x and y in either integral. Also, we 
note that if x ^ 0, or x ^ oo, then, by the transformation y = xz: 

which by the transformation y — ► 1/y equals 

1 r dy 1 
xj (y + lf^/y 

which can now be returned to 



(B.4) 



■"OO 



dy fx 



Jo (x + y) 2 \ y 
Hence if x ^ then 

f°° dy fy__ f 00 dy fx_ 
Jo (x + y) 2 \x Jo {x + yfyy' 
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(B.5) 



(B.6) 



However, if we break the integration of y from 



roo re poo 

/ to / + / , (B.7) 

JO J0 Je 



we see that 



/ dx I 

Jo Jo 



- dy hj r j r d y fx 



dx 



{x + y) 2 \l x J ' J ( x + y ) 2 \J y 



Jo dX Io {x + yfi}fx {~y) + lo {x + yfi^x {y) 

(B.8) 

The first term is zero by relabelling x y in one of the integrals, and 
the second term, by writing x — > ex, and y — > ey 



\2 



= a constant c > (B.9) 

since y > x throughout the entire range. 

So we have proved (using a similar proof at x = oo) that 



L 



dy i!J r (h! ^ = a6(x)-a6(x-oo). (B.10) 



'o (x + y) V x Jo (x + y) V V 
This model can be rigorized in the manner of 8 function theory, but 

since we offer it simply to clarify the situation for the reader we will not 

give the rigorous mathematics. The model is a faithful representation 
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of what goes on in equation (|4.2|) , in this particular input, if short of 
the algebraic complexity of the 3D Navier Stokes equation. This means 
that if we assume Kolmogoroff is true and substitute k~ T in equation 
( |1.15|) we can expect the non- linear terms to cancel at T = 11/3. 
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c 

In this appendix we discuss in detail the long time decay of $k(t), that 
follows from equation (|XT3|) To orient oneself, consider what the terms 
in lap look like if the solution were an exponential decay. The two 
"Boltzmann" terms decay like: 

$ k (t) and$ H (t)$j(t) (C.l) 

i.e. e^ kt and e -fa-i+"i)\ 
Further, think of what would happen if cjk were just viscosity, i.e. k 2 . 

Then the locus of j when 

^ 2 = (k-j) 2 + j 2 (C.2) 

is a d-dimensional sphere with k as its diameter, from Pythagoras's 
theorem. For j inside the sphere 

and for j outside the sphere, 

< ^k-j + Wj- (C.4) 
It becomes now very clear that the contribution of j's within the sphere 
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decays slower than $k(£) and cannot, therefore, solve eq. (|3.l3| ) for 
long times. The above discussion is limited to decay rates that are 
proportional to k 2 . It is well known that in non-trivial cases the decay 
rates are proportional to k^, with \i ^ 2, and in most cases (turbulence 
being an exception) fi > 1. It is easy to verify that if fi > 1 the 
minimum of cjk-j + ujj occurs at j = k/2 and the value of cjk-j + cjj at 
the minimum is 2 1 ~ M u;k- The only difference from the case of fi = 2 
is the shape of the region where tOk-j + tOj < tOk in which the vector 
k plays the role of an axis of symmetry and its end points are on 
the boundary of the region. The conclusion is that the single mode 
decay that can be a quite reasonable description of $k(^) for relatively 
short times (cu^t < 1) cannot describe the long time decay. Since the 
argument depends on whether fi is larger or smaller than 1, we will 
discuss in the following the two cases separately. 

a. The situation /i > 1 is generic, it appears in KPZ, dynamics at 
the transition of the 4 theory etc. We expect that, for long times, 

00 

$k(*) = X>(AO/»(fc*t) (C5) 

i=l 
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such that 



J£w =0 (c ' 6) 



and the /i(0)'s are chosen to be 1. If the series in ( |C.5| ) converges for 
t = 0, then 

00 

$>*(*) = (C7) 

but there is no reason to assume a scaling function of the form <j) q f(u q t). 
In fact, since the cases we consider have D\& ^ for k ^ 0, this is 
actually impossible, because of the exact relation ( |5.lY| ) that is also 
obeyed by our approximation 

lim « = - M . (C.8) 

t^o+ at 

This is possible with the single scaling function form only if cu^ = 
jTnh^ - This is not normally the case, cjk does not scale as ko/0k- Our 



aim is to obtain the functional form of fi(x). We assume first that the 
non-linear term is not relevant at large t, and arrive at a contradiction. 
If it is not relevant, eq. (|3."13| ) gives 



dt 



+ |z/ k + J A<%fc)0jj$ k = O (C.9) 
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and the solution is 



Mt) = 0ke- Wk * (CIO) 



with 



^ k = ^k + J A (1) (j, k)<f> } dj (C.ll) 

but since the cjk obtained here behaves like /c M for small fi, with /i > 1, 
such a solution is not consistent with the assumption that the non- 
linear term is not relevant. By the same reasoning, leading to the 
conclusion the e _a;k * is not consistent as the description of the long time 
decay, it becomes clear that any decay faster than an exponential for 
fix) is not acceptable. For a decay that is slower than an exponential 
the derivative can be neglected compared to the function, for large x 
and eq.( |3.13|) may be approximated for large co^t by 



+ J A%\Wjdjj$ k (t)- J A%\A0$j(t)$k-j(*)=O. 

(C.12; 
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Using the long time leading order behaviour we find 

+ J A (1) (j^)0jdjjtfViM) 

- J A^( J: ^^V^j/iCwj^/iC^-jt) = 0. (C.13) 

Consider first the one dimensional case. We assume, without loss of 
generality, that k is positive and break up the integral on the left hand 
side of eq.( |C.l5|) into two contributions, 



A (2 VjMVi(^)/i(^-j*) dj = 

A (2) (j, fc)^ ( M-j/i(^*)/i(^-J*) d J 
+ fA^^j^^^fi^hiu^dj, (C.14) 

J\0M 



f 

'o 



where [0,k] is the set complimentary to the segment [0,k]. Because 
of the form of cjk, it is obvious that the second contribution on the 
right hand side of ( |C.14| ) decays faster than the first contribution. The 
leading order decay in eq.( |C.13|) will be thus obtained by equating 

= jf * A< 2 >0, k)^ 1] ^fi(^t)fi(^t) dj (C.15) 
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Since this equation must hold for all long enough times it is obvious 
that the only way to satisfy the equation is to have /(u^t) oc e -7 ^ 1 
where 7 is a numerical constant. The reason is that within the range 
< j < k 

|k| = |j| + |k-j| (C.16) 

so that a single time decay can be brought out of the j integral and 
made equal to the decay of &k(t). (In fact, the situation is a little 
more subtle than the above description, because regardless how large 
is the time, within the range of integration there are always j's and 
(and k — j's) such that u^t < 1 and there the long time form for the 
two point function is unjustified. This region of small j's decreases 
however with t and in the generic case where there are no convergence 
difficulties at j = this becomes unimportant). The coefficients ip^ 
obey the equation 

|^ + J Wj}^"^ A {2 \j,k)^\ = 0. (C.17) 

Comparing equation ( |C.lY| ) with the static equation for 0k (equation 
3.1 til ), or the approximate form in the inertial range (equation 0) that 
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yields the leading order behaviour that reads in terms of and A^ 2 ^ 

poo poo 

/ A%\ k)^ k - / A%", A:)0j0 k _j = 0, (C.18) 

J oo J oo 

we find that to leading order in k 

= C0 k , (C.19) 

where 

C= ■f-°° A() ^ k -J dj ' (c.20) 
J Q k A^(j,k)cf>^dj 

can be shown to be independent of k. The faster decaying contribu- 
tions to the A^ 2 ) integral, coming from j's outside the segment (0, k) 
will be matched against the faster decays present in <&k(t) (eq. |U.5|) . 
In more than one dimension the situation is a bit more subtle. The 
reason is that the long time decay of &k{t), exp(— ^kt 1 ^), can be 
matched against the slowest decay in the A^ 2 ^ integral exactly as in 
the one dimensional case but that slowest decay is contributed by a 
set of j's that is of measure zero, i.e. j = ak with < a < 1. To 
obtain the correct behaviour we need to consider not only the set of j's 
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where the decays can be matched exactly but also the nearby vicinity 

|j| + |k - j| - k < B- 1//i ( 7 t 1/ ^)- 1 = e, (C.21) 

where e is small compared to k (this is the basic long time condition 
Co>k£ >> 1), such that the difference in the decays is quite small. (The 
constant B is given by = Bk^). It is interesting to note that the 
effective region of integration described by eq. (|C.^1|) is the interior 
of an ellipsoid of revolution 

(J '-f )2 + f = 1, (C.22) 
a z b z 

where j'ii is the component of j in the direction of k, and j± is the part 
of j perpendicular to k, 



a = £+* and b = J ^^ . (C.23) 



Equation (|C.15|) will be solved now by writing 

$ k (t) = CMuktfe-^ 1/fl : (C.24) 

where C is a dimensionless constant and restricting the integration 
to the effective region described above. We find that the form given 
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above for $k(A) is adequate to leading order provided [£2 



(a25) 

The constant C can be obtained explicitly as in the one dimensional 
ratio of two dimensionless integrals but it is not of much 
importance. 

b. The treatment presented above relies heavily on the assumption 
that fi > 1. This is generic yet the most important, perhaps, in the 
class of systems we study here is the noise driven Navier- Stokes, which 
has fi = 2/3. To treat that particular system, we go back to equation 
( p7T3|) and write 

$kW = k e- FkW (C.26) 



to yield 

dr k 

dt 



-Vk - f A (1 Vj d 3 j + 

f A (2) ^k z j e -r J (t)-r k _ J (t)+r k (t) d 3^- (c.27) 
J &k 

We subtract the above from the static equation ( |Q| ) (divided by 0k) 
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and taking into account that ko is zero for finite k, we obtain 

d£k = f A (2)0j0k-j A _ e -r j( t)-r k _ jW+ r kW A ^ (C 28) 

dt J k V / 

We assume now there is a scaling solution, r k (£) = V^u^t), and since 
cjk and 0k are each a power of k, we can scale out t by the transfor- 
mation k — > (BtY'^k, whereupon we find 

3 fc i /3 ar(fc) = i / A (2) (kJ)j -n/3 / jk-j| \ - 11 / 3 



X 



2 dk B 2 J k 2 / 3 V k 

(l - e -r(i)-r(|k-j|)+r(k)^ ^ (C 2g) 

Note that ^/i? 2 is dimensionless. 

We consider now (|C.29|) in the regime of very large /c, and assume 



that r(/c) oc k v ', with perhaps some logarithmic correction. 

We now split the integral on the right hand side of ( |C.29|) into con- 



tributions coming from regions where j (or |k — j|) is small compared 
to 1 and the region where both are large compared to 1. Consider 
first the region where both j and |k — j| are large. In that region, the 



term in the exponent on the right hand side of ( |C.29|) is proportional 



to —\j v + |k — }\ v — k^]. The discussion of the first part of this section 
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applies and the conclusion is that 

v < 1. (C.30) 

Consider next the region where j is small compared to 1. Here we 
cannot say that is proportional to j v ' . In fact, it may start off 
as j a for small j, and cross over to f for j's of the order of 1. In 
any case, consideration of the small j dependence of the integrand in 
equation ( |(J.29|) reveals that it is enough to have a > 2/3 in order to 



remove the convergence difficulties that may be caused by the j 11//3 
factor. Now for j's that are less that 1 we can approximate 

r(|k-j|) - r(k) = ^c^ + D^ + E^^^k-, (c.3i) 

where the constants C, D and E are obtained from the expansion. 
We will assume now that v is actually less that 1 in order to study 
the consequences of that assumption. In the case that v < 1, the 
approximation presented by eq.( |C.5l|) is indeed justified, because the 
term on the right hand side is small due to the fact that k >> 1. Note 
that this is not true for v = 1 where for j of order 1 we are left on the 
right hand side of ( |(J.31| ) with an expression that is of order 1. 
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The conclusion of all this is that if v < 1 then for large k 

r(j) + r(|k-j)-r(k)>o, (c.32) 

apart, perhaps, from a very small region of j's near the origin depend- 
ing on the actual value of a and is smaller than , i.e. vanishingly 
small as k becomes large, and therefore of no consequence to the in- 
tegral. The inequality holds for large j and large |k — j| trivially, 
because v < 1 and for j smaller than 1, the left hand side of (|C732 



may, in fact, be approximated by r(j). An approximation that as 
stated above can fail, depending on the value of a, only for very small 
j's such that j < k~ e , with 6 > 0. The conclusion now is that the 
integrand in eq.( |C.^91) is positive definite and therefore contributions 



from different regions sum up with no cancellation, (note that it is 
positive definite only for large k and v < 1). The leading order k 
contribution from the small j integration, which we denote by /<(/c), 
can be easily shown to scale like 

I<0) oc /c 2/3 . (C.33) 

This follows since for small j, A^ 2 ^(k, j) is proportional to Because 
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contributions are summed up, it is clear that the leading k dependence, 
I(k), obtained from the full integration must be such that 

I{k)>I < (k). ( C - 34 ) 

Consequently, since I(k) oc /c^ 2 / 3 , as can be easily seen from the left 
side of flU.^M ) we arrive at 

^ > 4/3. (C.35) 

that is a contradiction to our assumption that v < 1. Thus the as- 
sumption v < 1 is not consistent and we are left with v = 1 as the only 
possibility. Hence we have shown that, to within a possible logarithmic 
correction 

T(x) oc |4 (C.36) 
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D 



We have noted that the expansion parameter of the series is —jrk and 
in order for the series to scale, to is determined by this relation. Thus 
in KPZ at any order of accuracy, the equation governing the index is 
Z(T) = and Z is calculated using an to whose power is determined 
by 

M 2 <t)k d uj- 1 = 0(1). (D.l) 

The front factor B in = Bk^ does not enter the Z(T) of ( |b.3|) 
and this is crucial for the self consistency of the theory. Of course we 
have calculated numerically using ^(T) and get, say, T2. When we go 
as in flO ) to 



Z 2 + Z 4 = 0, (D.2) 

the value of T will certainly change for KPZ (but apparently not for 
KNS, for KNS is deficient in dimensions). One can only check that 
Z4 makes little difference by actually calculating it, which is not at all 
easy from sheer algebraic size, though we repeat the comment earlier 
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that the numerical value from agrees excellently with numerical 
experiment. 

It follows from this discussion that the only significance of an cu^ 
equation lies in the ability to calculate a front factor. The series is 
formally valid for any u, but we have one powerful constraint in the 
value of a scaling relation. We have chosen to use a time decay argu- 
ment to derive but other methods are present in the literature. 



We could study the Hermite operator and looked at, in harmonic os- 
cillator language, the first excited state above Po i.e. HPq and try to 
fit that [123]. However, all the equations used have Galilean invariance, 



K /i k e lk vt (D.3) 

and the 'excited state' equations do not seem to have it, i.e. a denom- 
inator 

+ + uj-k-j — > + + u-k-j + iv • (k + j — k — j) (D.4) 

i.e. cjk + cjj + co-k-j is Galilean invariant; but cjj + U-k-j — <^k is 
not, and it is this combination that appears in the 'excited state' 
equations; they also do not always converge. A difficulty of the same 

61 



nature arises in ref |jj and if one uses the method of Herring [jZ4 
where the combinations cjj + uj-k-j appear in the determination of 
static properties, that are not affected by Galilean transformations. 
Another method is to argue that since the equations should be valid 
for any cu one can minimise the effect of variation of cu and differentiate 
( |1.15|) w.r.t. cjk- Yet another method is to maximise the "entropy" 
in the [h(r,t)} space [[T7j. All these remarks apply equally to Q and, 
whereas for cu^ one is only searching for a front factor equation, with 
£lkcj the whole complex plane presents itself. We feel that a major 
problem has been uncovered here, and although we have presented a 
practical solution to the problems, there is still a great deal yet to be 
uncovered. 
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Figure 2: The coefficient Z as a function of T for the case of turbulence. 
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